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We study the noise of the cotunneling current through one or several tunnel-coupled quantum 
dots in the Coulomb blockade regime. We consider the regime of weak (elastic and inelastic) 
cotunneling, and prove a non-equilibrium fluctuation-dissipation theorem which leads to a 
universal expression for the noise-to-current ratio (Fano factor). 



1 Introduction 

In recent years, there has been great interest in the shot noise in mesoscopic systems, because 
it contains additional information about correlations, which is not contained, e.g., in the linear 
response conductance. The shot noise is characterized by the Fano factor F = S/el, the dimen- 
sionless ratio of the zero-frequency noise power S to the average current /. While it assumes 
the Poissonian value F = 1 in the absence of correlations, it becomes suppressed or enhanced 
when correlations set in as e.g. imposed by the Pauli principle or due to interaction effects. 
In the present paper we study the shot noise of the cotunneling Hi current. We consider the 
transport through a quantum-dot system (QDS) in the Coulomb blockade (CB) regime, in which 
the quantization of charge on the QDS leads to a suppression of the sequential tunneling current 
except under certain resonant conditions. We consider the transport away from these resonances 
and study the next-order contribution to the current a (see Fig. 1). One might expect that the 
cotunneling, being a two-particle process, leads to strong correlations in the shot noise and to 
the deviation of the Fano factor from its Poissonian value F = 1. As it has been found recently, B 
this is indeed the case for the regime of strong cotunneling, i.e. when the cotunneling rate I/e is 
large compared to the intrinsic relaxation rate w m of the QDS to its equilibrium state due to the 
coupling to the environment, I/e S> W[ n . However we find here, that in the weak cotunneling 
regime, I/e w- m , the zero-frequency noise takes on its Poissonian value, as first obtained for a 
special case.H This result is generalized here, and we find a universal relation between noise and 
current for the QDS in the first nonvanishing order in the tunneling perturbation. Because of the 
universal character of this result (Eq. |l2|) we call it the nonequilibrium fluctuation-dissipation 
theorem (FDT) 13 in analogy with linear response theory. 

2 Model system 

In general, the QDS can contain several dots^ which can be coupled by tunnel junctions, the 
double dot (DD) being a particular example. Q The QDS is assumed to be weakly coupled to 
external metallic leads which are kept at equilibrium with their associated reservoirs at the 
chemical potentials m, I = 1,2, where the currents can be measured and the average current 
I through the QDS is defined by Eq. B Using a standard tunneling Hamiltonian approach, H we 
write 

H = H + V, H = H L + H s + H int , (1) 
H L=Y.Y. £ k4 k c tk , H s = £ p d ] p d p , (2) 

1=1,2 k V 



a The majority of papers on the noise of quantum, dots consider the sequential tunneling regime, where a 
classical description ("orthodox" theory) is applicable. u In this regime the noise is generally suppressed below its 
full Poissonian value F = 1. This suppression can be interpreted □ as being a result of the natural correlations 
imposed by charge conservation. 




Figure 1: The quantum dot system (QDS) is coupled 
to two external leads I = 1, 2 via tunneling barriers. The 
tunneling between the QDS and the leads is parametrized 
by the tunneling amplitudes Ti, where the lead and QDS 
quantum numbers k and p have been dropped for simplic- 
ity. The leads are at the chemical potentials fii,2, with an 
applied bias A/i = /ii — /i2- The eigenstates of the QDS 
with one added electron (N + 1 electrons in total) are in- 
dicated by their energies Ex, E% , In the cotunneling 

regime there is a finite energy cost A±(l,N) > for the 
electron tunneling from the Fermi level of the lead I to 
the QDS (+) and vice versa (— ), so that only processes of 
second order in V (visualized by two arrows) are allowed. 



v=y,(Pi + D \), a = E T ik P 4k d P . ( 3 ) 

1=1,2 k,p 

where the terms Hi and H$ describe the leads and QDS, respectively (with k and p from a 
complete set of quantum numbers), and tunneling between leads and QDS is described by the 
perturbation V. The interaction term H mt does not need to be specified for our proof of the 
universality of noise. The iV-electron QDS is in the cotunneling regime where there is a finite 
energy cost A±(l,N) > for the electron tunneling from the Fermi level of the lead I to the 
QDS (+) and vice versa (— ), so that only processes of second order in V are allowed. 

To describe the transport through the QDS we apply standard methods! and adiabatically 
switch on the perturbation V in the distant past, t = to — > — oo. The perturbed state of the 
system is described by the time-dependent density matrix p(t) = e -iH(t-to) p e lH ^~ to \ with po 
being the grand canonical density matrix of the unperturbed system, po = Z~ 1 e~ K ' kBT , where 
we set K = Hq — J2i^i^i- Because of tunneling the total number of electrons in each lead 
Ni = J2k c \k c ik is no longer conserved. For the outgoing currents Ii = eNi we have 

Ii = ei [V, Ni] = ei(D\ -£>/)• (4) 

The observables of interest are the average current I = I2 = —I\ through the QDS, and the 
spectral density of the noise Swiu) = J dtSw (i) exp(iut) , 

Ii = Trp(0)4 S iv (t) = ReTrp(0)^(t)^(0) , (5) 

where 5Ii = ij — ij. Below we will use the interaction representation where Eq. || can be rewritten 
by replacing p(0) -> p and I t (t) -> U*(t)Ii(t)U(t), with 



U(t) = Texp 



-1 



dt' V(t') 



00 



(6) 



In this representation, the time dependence of all operators is governed by the unperturbed 
Hamiltonian Hq. 



3 Non-equilibrium fluctuation-dissipation theorem 

We note that the two currents // are not independent, because [Ii, I2] 7^ 0, and thus all correlators 
Sw are nontrivial. The charge accumulation on the QDS for a time of order A^ 1 leads to an 
additional contribution to the noise at finite frequency u>. Thus, we expect that for uj ~ A± 
the correlators Sw cannot be expressed through the steady-state current / only and thus I has 
to be complemented by some other dissipative counterparts, such as differential conductances 



Gui. On the other hand, at low enough frequency, co <C A±, the charge conservation on the 
QDS requires 5I S = {8I 2 + 5I±)/2 ~ 0. Below we concentrate on the limit of low frequency 
and neglect contributions of order of uj/ A± to the noise power. In the Appendix we prove that 
S ss ~ (u>/A±) 2 (see Eq. ^), and this allows us to redefine the current and the noise power 
as I = Id = (I 2 — h)/2 and S(ui) = Sdd(u). j] In addition we require that the QDS is in the 
cotunneling regime, i.e. the temperature is low enough, k^T <C A±, although the bias Ap is 
arbitrary as soon as the sequential tunneling to the dot is forbidden, A± > 0. In this limit the 
current through a QDS arises due to the direct hopping of an electron from one lead to another 
(through a virtual state on the dot) with an amplitude which depends on the energy cost A± of 
a virtual state. Although this process can change the state of the QDS (inelastic cotunneling), 
the fast energy relaxation in the weak cotunneling regime, Wi n 3> I/e, immediately returns it to 
the equilibrium state (for the opposite case, see Ref.H). This allows us to apply a perturbation 
expansion with respect to tunneling V and to keep only first nonvanishing contributions, which 
we do next. 

It is convenient to introduce the notation Di{t) = J_ dt' Di(t'). We notice that all relevant 
matrix elements, (N\Di(t)\N + 1) ~ e _iA +*, (N - l\Di(t)\N) ~ e iA ~ l , are fast oscillating 
functions of time. Thus, under the above conditions we can write Di(oo) = 0, and even more 
general, J+™dtDi(t)e ±lujt = (note that we have assumed earlier that u> <C A±). Using these 
equalities and the cyclic property of the trace we obtain the following results (for details of the 
derivation, see Appendix), 

00 

I = e J dt([Ai(t),A(0)]), A = D 2 d\ + d\d 2 , (7) 

—00 

00 

S(u>) = e 2 I dt cos(ujt){{A^(t),A(0)}) , (8) 



where we have dropped a small contribution of order u/A± and used the notation (...) = 
Trp (...). 

Next we apply the spectral decomposition to the correlators Eqs. and ||, a similar proce- 
dure to that which also leads to the equilibrium fluctuation-dissipation theorem. The crucial 
observation is that [Hq,N[] = 0, I = 1,2. Therefore, we are allowed to use for our spectral de- 
composition the basis |n) = \E n , Ni, N 2 ) of eigenstates of the operator K = Hq — J2i fJ>iNi, which 
also diagonalizes the grand-canonical density matrix po, p n = (njpo l n ) = Z^ 1 exp[— En/ksT]. 
We introduce the spectral function, 

A(u) = 2irJ2(pn + Pm)\(m\A\n)\ 2 5{LJ + E n - E m ) , (9) 

n,m 

and rewrite Eqs. [7| and || in the matrix form in the basis |n) taking into account that the operator 
A, which plays the role of the effective cotunneling amplitude, creates (annihilates) an electron 
in the lead 2 (1) (see Eqs. || and 0). We obtain following expressions 



I(Afi) = etanh 



A/i 



2 2k B T} A ^' (10) 
S(w,A/i) = i-Vj4(A (1 ±ai). (11) 



b We note that charge fluctuations, 8Q(t) = 2j_ dt'5I s (t'), on a QDS are also relevant for device applications 
such as SET. EH While we focus on current fluctuations in the present paper, we mention here that in the cotunneling 
regime the noise power {SQ 2 )^ does not vanish at zero frequency, (8Q 2 )u=o = 4uj~ 2 S ss {l>j)\u^o 7^ 0. Our formalism 
is also suitable for studying such charge fluctuations; this will be addressed elsewhere. 



We note that because of additional integration over time t in the amplitude A (see Eq. |7|), 
the spectral density A depends on fj,% and \i2 separately. However, away from the resonances, 
u) <C A±, only A/i-dependence is essential, and thus A can be regarded as being one-parameter 
function, f] Comparing Eqs. 10 and |TI], we obtain 

" A/u ± uj~ 



S(uj,A/j) = ^^coth 
z ± 



2k B T 



I(Afjt±w) (12) 



up to small terms on the order of u/A±. This equation represents our nonequilibrium FDT for 
the transport through a QDS in the weak cotunneling regime. A special case with T, uj = 0, 
giving S = el, has been derived earlier.! To conclude this section we would like to list again the 



conditions used in the derivation. The universality of noise to current relation Eq. |12| proven 
here is valid in the regime in which it is sufficient to keep the first nonvanishing order in the 
tunneling V which contributes to transport and noise. This means that the QDS is in the weak 
cotunneling regime with u, k B T <C A±, and I/e <C w m . 
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Appendix 

In this Appendix we present the derivation of Eqs. ^ and ^. In order to simplify the intermediate 
steps, we use the notation 0(t) = f^^dt'O^') for any operator O, and O(0) = O. We notice 

that, if an operator O is a linear function of operators D\ and Dj, then O(oo) = (see the 
discussion in Sec. ||). Next, the currents can be represented as the difference and the sum of I\ 
and I2, 

Id = (I 2 -h)/2 = ie(Xi-X)/2, (13) 
Is = (/1 + / 2 )/2 = ie(y f -Y)/2, (14) 

where X = D2 + D\, and Y = D\ + D2. While for the perturbation we have 

v = x + x f = y + y f . (15) 

First we concentrate on the derivation of Eq. |7| and redefine the average current Eq. || as / = Id 
(which gives the same result anyway, because the average number of electrons on the QDS does 
not change I s = 0). 

To proceed with our derivation, we make use of Eq. || and expand the current up to fourth 
order in T lkp : 

t 
I = i J dt J dt'(i d V(t)V(t')V{t')} -i J dt(VI d V(t)V(t)) + c.c. (16) 

—00 —00 —00 

Next, we use the cyclic property of trace to shift the time dependence to Id- Then we complete 
the integral over time t and use Jrf(oo) = 0. This procedure allows us to combine first and second 



term in Eq. 16 







-i J dt([I d V + VI d ]V(t)V(t)} + c.c. (17) 



c To be more precise, we neglect small tj-shift of the energy denominators A±, which is equivalent to neglecting 
small terms of order oj/A± in Eq. |ll|. 



Now, using Eqs. [l3| and [15] we replace operators in Eq. |l7] with X and in two steps: I = 
e/_ OG dt([X+Xt - XX]V{t)V{t)) + c.c., where some terms cancel exactly. Then we work with 

V(t)V(t) and notice that some terms cancel, because they are linear in cik and cL. Thus we 
obtain I = e dt{[X* X^ - XX][X^(t)X^(t) + X(t)X(t)]) + c.c. Two terms XXXX and 
X^XTX^XT describe tunneling of two electrons from the same lead, and therefore they do not 
contribute to the normal current. We then combine all other terms to extend the integral to 

+00, 

00 

I = e J dt(X\t)X\t)XX - XXX\t)X\t)) (18) 

—00 

Finally, we use dtX (t)X (t) = - /f^ dtX(t)X(t) (since X(oo) = 0) to get Eq. @ with 
A = XX. Here, again, we drop terms -Of-D} and D2D2 responsible for tunneling of two electrons 
from the same lead, and obtain A as in Eq. [7]. 

Next, we derive Eq. ^ for the noise power. At small frequencies u <C A± fluctuations of I s are 
suppressed because of charge conservation (see below), and we can replace I2 in the correlator 
Eq. H with Id- We expand S(oj) up to fourth order in T/fc p , use f_™ dt Id{t)e ±%U)t = 0, and repeat 
the steps leading to Eq. 17. Doing this we obtain, 

00 

S(u) = - [ dtco S (Lot)([V(t),i d (t)}[Vj d }). (19) 



Then, we replace V and l d with X and X*. We again keep only terms relevant for cotunneling, 
and in addition we neglect terms of order u/A± (applying same arguments as before, see Eq. [20|) . 
We then arrive at Eq. ^| with the operator A given by Eq. |7| 



Finally, in order to show that fluctuations of I s are suppressed, we replace Id in Eq. |19| with 
I s , and then use the operators Y and Y ' instead of X and X^ . In contrast to Eq. ^ terms such as 
ytytyy do not contribute, because they contain integrals of the form /f^, dt cos(ujt)Di (t)Di> (t) = 
0. The only nonzero contribution can be written as 



e 2 oj 2 



S ss (lo) = — - \ dtco S (ut)([Y\t),Y(tW\Y]} , (20) 



where we have used integration by parts and the property Y(oo) = 0. Compared to Eq. || this 
expression contains an additional integration over t, and thereby it is of order (lu/ A±) 2 . 
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